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A pair of conjugate observables, such as the quadrature amplitudes of harmonic motion, have 
fundamental fluctuations which are bound by the Heisenberg uncertainty relation. However, in a 
squeezed quantum state, fluctuations of a quantity can be reduced below the standard quantum 
limit, at the cost of increased fluctuations of the conjugate variable. Here we prepare a nearly 
macroscopic moving body, realized as a micromechanical resonator, in a squeezed quantum state. 
We obtain squeezing of one quadrature amplitude 1.1 ± 0.4 dB below the standard quantum limit, 
thus achieving a long-standing goal of obtaining motional squeezing in a macroscopic object. 


The motion x(t ) = X\ (t) cos(cu m t) + X 2 (t) sin(cj m t) of 
a harmonic oscillator having the natural oscillating fre¬ 
quency co m can be described by the quadrature ampli¬ 
tudes X\ and X 2 which have slow fluctuations. The fluc¬ 
tuations, presented in units of the quantum zero-point 
fluctuation amplitude x zp , satisfy the Heisenberg uncer¬ 
tainty relation XX\XX 2 > 1. One of the two can be 
prepared (=squeezed) below the value 1, at the expense 
of increased fluctuations in the other quadrature. In op¬ 
tics, squeezing of laser light was observed in early 80’s 
m, not long after the possibility was realized. 

It has been a formidable challenge to obtain squeez¬ 
ing in the motional state of a macroscopic object. The 
possibility of squeezing in the oscillations of massive grav¬ 
itational antennae has been hypothesized a long time 
ago om, but technological limitations are too severe 
for experimental realization. Other motional quantum- 
mechanical phenomena, on the other hand, have recently 
been experimentally demonstrated mi in micromechan¬ 
ical resonators. The latter systems are nearly macro¬ 
scopic in physical size, and therefore they provide an ideal 
test system for treating the borderline between quantum 
and classical. Of particular interest for these studies has 
been the cavity optomechanics setting coupling electro¬ 
magnetic cavity mode and the oscillator motion [7] . Out¬ 
put of squeezed light IMEQ] was recently observed, but 
this does not yet imply the oscillator mode is squeezed. 

Here we report the first realization of squeezing of the 
motional state of a nearly macroscopic body, realized as 
a micromechanical resonator measuring 15 microns in di¬ 
ameter. We utilize the novel idea of dissipative squeezing 
mm , see Fig. 0a, where the system is allowed to cool 
towards a squeezed low-energy state. This method has 
the great advantage of being able to create unconditional 
squeezing in the steady-state. This contrasts to many 
other plausible methods of squeezing generation [[T4H2Q] . 
Our approach is closely related to the quantum nonde¬ 
molition measurements EIH22] which, however, are not 
able to generate true squeezing without feedback. At 
this point we mention that classical squeezing of thermal 
noise is routinely observed in mechanical systems [241126] . 

The mechanical element (Fig. 0b) is a drum resonator 


basically similar to Refs. [27j |28| The 120 nm thick 
A1 drum fabricated on a quartz substrate is connected 
through a narrow 70 nm vacuum gap to one end of the 
cavity. The latter consists (Fig. |T] c) of a meandering 
superconducting A1 strip. This kind of microwave op¬ 
tomechanical systems are conveniently described using a 
model involving lumped electromagnetic elements. The 
interaction between the microwave cavity and the me¬ 
chanical vibrations is given by the energy gonpx , where 
np is the number of photons externally applied by mi¬ 
crowave pump(s), and go is the radiation-pressure cou¬ 
pling energy. Because usually go is much smaller than 
other energy scales, a large np 1 is needed to effec¬ 
tively enhance the interaction up to a value G = goftnp- 
If the pump is applied at the frequency = uj c — uj m 
(the red sideband), the physics leads to sideband cool¬ 
ing of the mechanical vibrations, possibly down to the 
quantum ground state [28], 29- 

The scheme of Ref. m requires two pump microwaves, 
one applied at the red sideband, and the other at the blue 
sideband frequency uj+ = uj c + (Fig. |T] d). The two 

pertinent effective couplings are called G_ and G+, re¬ 
spectively. This setup can be described as sideband cool¬ 
ing of a Bogoliubov (BG) mode [TTJ 30] , which in labora¬ 
tory frame corresponds to cooling the mechanical mode 
towards a squeezed vacuum state. The BG mode is de¬ 
fined by the annihilation operator [3 = b cosh r + sinhr, 
obtained from the creation and annihilation operators b t, 
b of the mechanical resonator. An arbitrary squeezing ra¬ 
tio r can be selected by tuning the ratio of the two pumps: 
tanhr = G + /G_. The BG mode - cavity system is de¬ 
scribed by the Hamiltonian H = Gp {ft f3 + aft) with 

the coupling energy Gp = yjG 2 _ — G\. Here, a\ a are 
the creation and annihilation operators of the cavity. Al¬ 
though a value G_/G+ > 1 would give rise to a high 
squeezing ratio, the effective coupling needed for cooling 
this mode would vanish. Hence, an optimum value is 
typically around G_/G + ^1.5. 

The measurements are carried out inside a Bluefors dry 
dilution refrigerator in the temperature interval 7 mK ... 
200 mK. The cavity is first characterized having the fre- 
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FIG. 1: Setup of the microwave optomechanical squeezing experiment, (a), Idea of dissipative squeezing. The initial fluctuation 
amplitudes of the quadratures are marked with dashed lines, and the final ones with solid lines. The green circles denote the 
quantum ground state. In sideband cooling (left), the initial fluctuations uniformly cool towards the ground state. Right, in 
a suitably engineered system, cooling can be quadrature-dependent, hence leading to one quadrature becoming squeezed, (b), 
(c), Optical micrographs of the micromechanical device and of the cavity. There are two drum resonators connected to cavity, 
however, only one of them (within the dashed rectangle) is operated. The cavity is asymmetrically coupled to the input port 
(left), and to the output port (right) for transmission measurements, (d), The frequencies involved in the scheme. The cavity 
is pumped by two nonequal-amplitude microwave tones at the sideband frequencies u± = uj c ±cj m . (e) Example of the thermal 
motion signal measured at the refrigerator base temperature, with the coupling G~/2tt ~ 1.3 kHz. 


quency cj c /27t ~ 6.9004 GHz, and coupled to the trans¬ 
mission lines with the decay rates kei/ 27t ~ 50 ± 5 kHz, 
keo/ 27t ~ 270 zb 30 kHz through the input and output 
ports, respectively. The internal losses are characterized 
by the rate ki/2tt ~ 330 =b 40 kHz, and the total cavity 
damping rate is ft = kei ± keo + ft/ — (27t) -650± 10 kHz. 
We operate in the good cavity limit cj m /ft ~ 20 1, a 

prerequisite for efficient sideband cooling and squeezing 
generation. 

The mechanical resonator is first characterized using a 
single pump tone at the red sideband. We choose very 
low pump powers such that the cavity back-action damp¬ 
ing rate = 4 G 2 _/e is much smaller than the intrinsic 
linewidth of the mechanics. The emission at the cav¬ 
ity frequency then shows the usual thermal motion peak 
at the motional sideband at a frequency uj m /2 tt ~ 13.032 
MHz above the pump. We obtain 7 m /27r ~ 330 Hz cor¬ 
responding to the Q-value Q m ~ 3.9 x 10 4 from the data 
as in Fig. [l] e (y_ ^ 12 Hz was subtracted from the fit 
result). 

An important benchmarking for cavity optomechanical 
experiments is how well the mechanical mode thermalizes 
to the temperature T of the refrigerator. Here we observe 
the linear temperature dependence expected (as seen in 
Fig. 1 a) in equilibrium, ksT = down to ~ 25 

mK. Here, n ^ is the equilibrium thermal phonon number 
defined accordingly. In what follows, we operate at the 
minimum T where we know that the mechanical mode is 
at 25 mK corresponding to ~ 40 at low pump powers. 

We proceed with a series of further calibrations on the 
way towards demonstrating squeezing. Next we perform 
a regular sideband cooling experiment (G+ = 0 ), see 
e.g. Refs. [3lti35] . For calibrating G versus generator 
power, we study the peak width + 7 - as a function 


of power at modest values of G_ ft. The most critical 
step, which also will account for most of the final im¬ 
precision, is to calibrate the system gain at the detection 
side. Each sideband cooling spectrum (about 20 curves at 
different power) are simultaneously fitted to theory with 
the same gain, using the G_ and n ^ just calibrated. For 
details, see the Supplementary [30] . We show examples 
of the sideband cooling spectra by the black symbols in 
Fig. [2]b-e, overlaid with theoretical predictions from the 
standard formalism using input-output theory [30] . For 
the plot, we have subtracted a large background level due 
to the amplifier noise, hence displaying only the signal 
part due to the sample. We also find that the mechani¬ 
cal mode cools down to a thermal occupation n m ~ 0.38 
(n m = 0 corresponds here to the ground state). The 
double peak seen in Fig. [2] e signifies the onset of the 
strong-coupling regime when G_ ^ ft. The data is plot¬ 
ted in dimensionless units (quanta) which is the natural 
unit from the theory point of view (W/Hz units are ob¬ 
tained by multiplying by Huj c ). 

Given that we can cool the drum motion close to the 
ground state provides a promising starting point for cre¬ 
ating squeezed motional states. We switch on the blue 
pump while keeping the red on, obeying G_ > G+ 
for stability. This creates a certain BG mode depend¬ 
ing on the ratio G_/G+. In order to ascertain which 
BG mode we have, we calibrate the input line attenu¬ 
ations separately for both pump generators. We select 
7 ± = 4G±/ft < 7 m and adjust the powers such that 
we obtain equal response due to either pump. The im¬ 
precision is estimated to be ±0.2 dB which is also the 
imprecision for constructing a given BG mode. 

Next we discuss a specific BG mode obeying G_/G+ ~ 
1.52 which is expected to represent a choice close to opti- 
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FIG. 2: Cooling the Bogoliubov mode. In all the panels, black 
color refers to regular sideband cooling (i.e., G+ = 0) used as 
calibration, whereas green refers to the BG mode experiment, 
(a), Thermalization in equilibrium. The left scale gives n^, 
as well as the area of both the sideband cooling peaks, and 
of the BG mode peaks (both in a.u.), whereas the right hand 
side scale is n^ G . (b)-(e), Output spectrum during sideband 
cooling (black), or under BG mode conditions (green). The 
solid lines are theoretical curves. The n m values refer to the 
G+ — 0 case. 


mum. At the lowest pump powers, the back-action cool¬ 
ing is negligible and we reveal the bare BG mode un¬ 
dressed from the cavity. Under this condition, the equi¬ 
librium BG mode occupancy n^ G is expected to follow a 
linear temperature dependence same way as the mechan¬ 
ical mode, but with an elevated temperature m- More 
relevant than however, is the agreement of the spec¬ 
tra with theoretical prediction, which here is connected 
to the area under the Lorentzian BG mode peak. We test 
this in Fig. [2] a and observe an excellent agreement to the 
theory. The green solid line is an expectation based on 
the calibrated G_/G + ratio, and on the G + = 0 data. 

Next, we increase the effective couplings. In order to 
mitigate possible gain drifts, we repeat a sequence of 
short measurements of the amplifier noise background, 


sideband cooling, and BG data. Plotted in conjunction 
with the sideband cooling data, in Fig. [2] b-e we display 
the BG mode spectra in green. The theoretical plots 
show a good agreement to the data. In Fig. [2] e, the 
BG mode curve is slightly shifted to the left probably 
because of an energy-dependent shift in the cavity fre¬ 
quency. When using fixed pump frequencies as here, the 
cavity can become slightly detuned at certain pump pow¬ 
ers. Here we have used as adjustable parameters the bath 
temperatures of both the cavity and the mechanics. Both 
baths heat up with the pump powers, for the BG mode 
up to ni ~ 0.74, ~ 86 (instead of we prefer 

to quantify the bath with n^) in Fig. [5]e, which is at¬ 
tributed to dielectric heating. These factors together set 
the limits for the cooling as well as squeezing. As one 
can anticipate, the baths heat up more if we apply both 
pumps instead of only one. 

In Fig. [2] b-e we label the fitted values of the BG mode 
occupation ubg which is a quantity analogous to n m in 
case of regular sideband cooling. Although ubg does not 
directly correspond to a physical temperature, a value 
Ubg < 1 as obtained in Fig.[2]e is indicative of squeezing. 
We can more precisely infer the amount of squeezing by 
evaluating the quadrature variances [30 . This way, we 
obtain that the mechanical mode is squeezed below the 
vacuum level, in Fig. [2] e by ~ 1.1 dB. 

The best verification of squeezing comes from the 
quadrature spectra which amount to tomography of the 
state. We digitally mix down the signal using the center 
frequency uj c of the pumps as a local oscillator (LO), 
hence making homodyne detection. The quadrature 
spectra show strong dependence on the LO phase. We 
identify the minimum emission as the X\ quadrature, 
and the maximum, offset by 90°, as X 2 . In Fig. [3] we 
show the corresponding quadrature spectra Si and S 2 
together with the total spectrum (the BG mode) and the 
cooling spectrum. We plot the quantities 2Si and 2 S 2 
for more conveniently presenting them together with the 
other two curves. The theoretical predictions overlaid on 
the data show an excellent agreement. In the best case 
of Fig. [3] d, the motion of the mechanical resonator has 
been squeezed about 1.1 dB below the Heisenberg limit. 
When varying the LO phase, we also observe an excellent 
agreement to theory (Fig. [4]). The bath temperatures are 
found to be slightly enhanced over the previous dataset 
in Fig. m mt We also find that this data agrees with 
a slightly shifted BG ratio G_/G + ~ 1.43 and which is 
attributed to drift, not directly measured, in the tunable 
filters at room temperature over the about one week after 
the calibration. 

We now discuss an issue which is critical for under¬ 
standing the quadrature spectra, namely, parametric ef¬ 
fects beyond the ideal optomechanical model. One can 
satisfactorily model the total BG mode spectrum (and 
get equal squeezing) without imposing any parametric 
modulation ED], however, the quadratures show much 
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FIG. 3: Squeezing inferred from the quadrature spectra. In all panels, blue and red refer to the cold and hot quadratures X\ 
and X 2 , respectively. Black and green refer to the regular sideband cooling, and the BG mode, respectively. The solid lines 
are theory curves. The pump powers are increased from (a) to (d) as marked in the panels while G-/G+ ~ 1.43 ratio is kept 
constant. The variances AXf are marked, and a value less than one implies squeezing below vacuum. The amplitudes of 
parametric modulation to the cavity are e c /2i: ~ 35,48,49,56 kHz from (a) to (d). 




G_/2tt (kHz) 

FIG. 4: Tomography and final results, (a), The quadrature 
data similar to Fig. 3 d, plotted at different LO phase values 
at 7r/10 steps from 0 to 7r/2, from bottom to top. (b), The 
Xi quadrature variance as a function of pump power. The 
blue region signifies squeezing below the quantum limit. The 
inset shows the quadrature variance from the data in (a). 


more phase dependence than predicted. In the scheme, 
mixing products of the pumps can appear both at 2uj rn 
as well as 2u; c , the frequencies most susceptible to cause 
parametric modulation via e.g. thermal effects or non- 
linearities [22j 30 . Modulation of uj c out of phase with 
respect to the pumps gives an excellent match to the 
quadratures with the values quoted in Fig. [3] In order 
to verify the existence of said parametric effect, we car¬ 
ried out a measurement where we substantially detuned 
both pumps from the sideband resonance, by ^ ±6ft, 
such that the center frequency stays at uo c [30]. This 
high detuning essentially eliminates optomechanical phe¬ 
nomena, but a possible field at 2uj c remains. This way 
we measure a spectrum consistent with a parametrically 
modulated oscillator with the correct phase. A possible 
explanation is nonlinear dissipation in the cavity [36], or 
a thermal effect. Although the parametric effects have 
a dramatic influence on the quadratures of the output 
spectrum, they only weakly affect the squeezing of the 
mechanics, in the present case we find a reduction by 
about 10 %. 

For the error analysis, we use a worst-case scenario of 
systematic errors from the calibrated parameters, and of 
direct statistical errors of the adjustable parameters. We 
find that the gain uncertainty is dominating. Because the 
shapes of the spectra are sensitive to most parameters, 
but squeezing is somewhat insensitive to any parameters 
in the present parameter region, we obtain a relatively 
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small imprecision of ^ 10 %. The final numbers are 
quoted in Fig. [4]b. 

In conclusion, we have achieved the reduction of the 
quantum fluctuations in a nearly macroscopic moving 
object below the ” sound of silence” level by 1.1 dB in 
one quadrature. Following the discoveries of ground-state 
cooling [5[ |28j ES] , observation of zero-point fluctuations 
in mechanical systems and entanglement [6], our 

work goes forward to establish the reality of another fun¬ 
damental quantum property in moving objects. 

Note added: while finalizing the paper, we became 
aware of a similar experiment [38^ . 
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EXPERIMENTAL SETUP 


We first mention how the microwave optomechanical device can be simply described using lumped electromagnetic 
elements (see Fig.j5|. The motion of the drum changes its capacitance C g (x). The total capacitance of the cavity can 
be summed up as a constant C and an x-dependent part C g (x). The frequency of this cavity hence is 


C y/HC+Cg{x))‘ 

Linearization gives the coupling 


9 o 


du c 


LJc_ 

2 C 



( 2 ) 


Modeling the device layout using electromagnetic simulation software gives the equivalent parameters C ~ 16 fF, 
L ~ 14 nH. It is beneficial to get as low C as possible in order to maximize coupling. Here the role of quartz substrate 
is critical because it has a low dielectric constant e r ~ 4 as compared to sapphire or silicon (e r ^ 10... 12). This 
difference amounts to a factor of two in C. 

A narrow vacuum gap ^ 70 ± 10 nm between the drum and a bottom electrode deduced from the device structure 
corresponds to dC g /dx ~ 120 nF/m, x zp ^ 5 fm, and hence allows for a relatively large single-photon coupling energy 
go/27r ~ 90 =b 20 Hz. We note that in the analysis, we are not directly using go, but the calibrations and analysis is 
using the derived quantities G±. 
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FIG. 5: Schematics of the microwave optomechanical device, (a), Model in terms of lumped electromagnetic elements, (b), 
Modeling scheme in the input-output theory. The cavity is exposed to noise from both the input (subscripts i) and output 
(subscripts 6) couplers, as well as from the internal losses (/). 


We use a pulse-tube powered Bluefors ’’dry” dilution refrigerator to cool the device down to about 7 mK. Inside 
the cryostat, the incoming pump signal is attenuated at all temperature stages by a total of about 40 dB (including 
cabling) in order to dampen the Johnson noise from higher temperatures. We used relatively small-valued attenuators 
at the two lowest temperature stages in order to avoid heating the refrigerator by the pump microwaves. However, 
one can check that the input noise n^, see Eq. (21), (considering an ideal setting) is only a vanishing amount above 
the vacuum level at ^ 7 GHz. The entire cabling is represented in Fig. [6] The signal from the sample is fed via 
mostly superconducting coaxial cables and two isolators into the low-noise amplifier at 4 Kelvin stage of the cryostat. 
A band pass microwave filter (BPF) is used to cut noise outside the isolator band. 

We will turn the discussion to the room-temperature instrumentation as shown in Fig. [7] At the input side, the two 
pump tones from two low-phase noise microwave sources are combined with a power splitter. A tunable notch filter 
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from Wainwright Instruments at the signal input at room temperature provides about 50 dB rejection ratio at the 
pump frequencies as compared to the center frequency cj c , and is meant to clean the phase noise of the pump sources. 
The phase noise would appear as an increased thermal noise ni to the cavity, limiting the performance. Without the 
filter, indeed, there is noticeable heating of the cavity at the highest used pump powers. 


I I room 
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FIG. 6: Setup for the cryogenic microwave measurements. 

The input signal is further split, and one branch is used to create the local oscillator (LO) at the center frequency 
which is phase-coherent with both the two pumps. The branch for LO creation is marked in Fig. [7] as the shaded box. 
The signal is squared which creates a tone at 2 cj c , then filtered around this frequency in order to block unintended 
tones at 2 cj c ± 2cj m . The frequency prescaler ADF5000BCPZ then creates the desired tone at uj c . It is attenuated to 
a suitable amplitude before being summed up to the signal from the sample, to be fed to the signal analyzer. 

At room temperature, the signal coming from the sample is first further amplified by 30 dB, then band-pass filtered 
within the bandwidth of 4 MHz around the cavity frequency. This filtering both acts as a antialiasing filter for the 
digitizer, as well as attenuates the pump tones by 15-20 dB to avoid saturation of the next stages. 

The vector network analyzer (VNA) is used to monitor the cavity response and hence enable for example correcting 
for occasionally observed drift in the cavity frequency. 

We use Anritsu MS2830A signal analyzer which digitizes the datastream at 2 MHz bandwidth for the duration 
of approximately 1 second, and automatically makes the conversion into IQ data which is sent to computer. At 
the computer, we first extract the phase of the LO marker signal from the time-domain data at about 10 ms time 
resolution, and use this information, first of all, to make small corrections to the data due to phase drift. We construct 
the quadrature spectra as follows. The IQ data is converted into real-valued time-domain data as 

I(t) = \IQ(t)\ cos(arg (IQ(t))), 

Q(t) = \IQ(t )| sin(arg(J Q(t))). 

These quadratures are Fourier-transformed. We then construct a set of phase-rotated frequency-domain datasets as 

Sv,e(w) = /©(cj) = cos (O)I(oj) — sin (O)Q(cj) . 

The final result of the processing is the quadrature spectrum |SV,©(^)| 2 - 
























The homodyne detection together with the fact that the peaks are centered at the LO frequency, cause that only 
one side of the quadrature spectra peaks are observable. For convenience, we however plot the data mirrored in about 
the center frequency. In the plots the data points do not appear symmetrically situated, however, which is due to a 
moving average filtering applied. We cannot know if the original quadrature spectra are symmetric, however, the fact 
that the total output spectrum (the BG mode) is symmetric about the center frequency, as well as the theoretical 
prediction, indicate they are symmetric. 

For obtaining the total spectra for the BG mode, or in case of sideband cooling, we directly process the spectrum 
as 


SVM = \St[IQ{t)]\ 2 . 


All the instruments are locked into a common 10 MHz Rubidium frequency standard. We noticed that the Ru clock 
improves the phase stability by about 20 dB as compared to using a quartz clock, substantially easing the digital 
processing. 



BPF o) c 


FIG. 7: Setup for the room-temperature microwave measurements. For discussion, see text. 


It is important to prevent the LO signal at uo c from entering the sample, because it would add a drive to the 
mechanics. We noticed that several loose connectors can be enough to allow for this tone to pass to the sample as a 
crosstalk through air. In this case we could see an enhanced temperature of the mechanics if the LO generation was 
switched on. In the final setup, this problem disappeared. 

Since we are not using a nearly quantum-limited Josephson parametric amplifier, the signal-to-noise is often much 
smaller than one. This leads to the issue that data recording takes long time, which causes the signal level to be 
sensitive to drifts of the total gain of the amplification chain. Because the vertical scale of the output spectra needs 
to be known accurately, it is imperative to make sure gain drifts are taken care of. To this end, while taking the 
quadrature spectra, instead of averaging one curve for a long time, we go through measuring a set of curves each with 
a different pump power. At a time, we record one curve for about 2 seconds. Before switching to the next one, we run 
two calibrations. The first is a measurement of the background, namely, switching all pumps off and recording the 
amplifier noise spectrum over the same band. This allows neutralizing gain drifts. The second is the regular sideband 
cooling which allows keeping track of possible drift in the input attenuation and thermalization. 
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THEORETICAL MODELING 


The cavity is described by the annihilation and creation operators a and a t, and the mechanical resonator similarly 
by b and We use the quadratures of the mechanics 


X x = + b 

X 2 = i (6+ - b) , 


(3) 


or the other way round: 


b=-(X 1 + iX 2 ) 
fet = i (X! - iX 2 ) . 

Involving an arbitrary quadrature angle 0, one can write 

X © = b t exp(i0) + 6exp(—iO). 

Similar relations of course hold for the cavity operators which we do not write explicitly. 


(4) 


(5) 


Basic Hamiltonian 

The Hamiltonian of the system consisting of the radiation-pressure coupled cavity and the mechanical resonator is 

H = cj c A^A + u m tfb + goA^Ax . (6) 

The uppercase cavity operators A and A^ involve both the quantum fluctuations described by the operators a and 
a as well as the rotating classical fields due to the two pumps: 

A = a- exp(— iw-t) + a+ exp(— iu+t) + a . (7) 

The quantities a+, a_ are the classical intra-cavity amplitudes. We can, without loss of generality, suppose that they 
are both real. A phase difference between the amplitudes would only shift the phase reference of the quadratures. We 
define the detunings of both pump frequencies from the cavity: 

A± = (8) 

In the experiment we operate very close to the exact sideband co-resonance condition A± = however, for the 

modelling, we relax this assumption and allow basically arbitrary A±. 

We transform to a rotating frame with 


U = ex.p(iuj' c ta^a) exp (iLj f m tb^b ). 


(9) 


At the sideband co-resonance, the arbitrary frequencies cJ c and c equal the cavity and mechanical frequencies, 
respectively. 

In accord to usual practice, we neglect higher-order terms in the fluctuation operators in what follows. This amounts 
to linearization owing to |a±| a. The coupling term becomes 


G16exp(i(— uo c + uj- — + G+6exp(i(— J c + cj + — 

+ G_6exp(i(— uj_ + uj' c — + G + 6exp(i(—+ J c — h 

+ G*_b ^ exp(i(— uj' c + + c + G^_b^ exp+ cj^)t)a-|- 
+ G-tf exp(i(— uj- + cj' c + c )cl [ + G+i^ exp(i(—+ uJ c + c jj f m )t)a ^. 


( 10 ) 


Here, the effective couplings are 


G± — go&zL 


(11) 
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We select ft — (cj + + u;_)/2, and = (cj + — )/2, simplifying Eq. ( |10| ) to 

G*_bft + G+6a +■ G_6a^ + G+&V + G+ba) exp(— 2i(jo' rn t)-\* 

+ G15aexp(— ‘liuoft) + G^^aexp(2icj^t) + G_&W exp(2icj^t). 


( 12 ) 


We are working in the good-cavity limit uJm/u ~ 20 >> 1, and hence we neglect the terms rotating at 2c ft. Relatively 
close to the sideband co-resonance, this becomes ~ 2cj m , and the rotating terms are suppressed by nearly 40 dB. We 
obtain that the Hamiltonian of the double-pumped system can be written as a beam-splitter between the cavity and 
the Bogoliubov modes: 


H = - 


A + + A 


aft ■ 


A: - A. 


UJm, 


-) tfb + a ] (G-b + G+rf) + a (G*_tf + G* + b) 


= Asa 1 "a + A A bft + Gp (a)ft + aft) , 
where the symmetric A s and asymmetric A^ detunings were defined, and 


Gs — \ G 2 _ — G\ 


(3 = b cosh r + sinh r 
tanhr = G + /G_ . 

The temperature of the BG mode is given as 

ubg = ftP = bftu 2 + (ft) 2 uv + b 2 uv + bbft 2 = + ft) 2 uv + 6 2 m; + (6^6 + l)(?x 2 — 1) = 

= 26^6 cosh 2 r — 6^6 + cosh r sinh r ((^) 2 + b 2 ) + sinh 2 r . 

At low effective couplings, the phase-sensitive correlations vanish, and we get 

ubg = b ft (cosh 2 r + sinh 2 r) + sinh 2 r . 


(13) 


(14) 


(15) 


(16) 


Parametric modulation 

On top of the ideal model in section [J we consider two additional effects, namely, a possible parametric modulation 
to the mechanics or to the cavity. We show below that the latter is critical for explaining the quadrature spectra, 
whereas the former does not play a significant role. The former was investigated in detail by Suh et al [23], who found 
that it can be a serious limitation in the double pump scheme if it is strong enough to cause parametric instability. 

In general, the double pump scheme is prone to issues due to unintended parametric modulation. Said briefly, this 
is because mixing products of the pumps at uj c ± uj m can appear both at 2cj m as well as 2cj c , the frequencies most 
susceptible to cause parametric effects. 


Mechanics 


In the scheme, there is an unavoidable parametric modulation resulting from the second-order optomechanical 

8 2 uj c _ ujn d 2 C a 


coupling with the coupling coefficient g 2 = 
The Hamiltonian Eq. © is added with 


2 c dx 2 


H 2 = ^g 2 A ] Ax 2 = hg 2 A^A ( b 2 + 2 b ] b + (& f ) 2 ) . 

By the ansatz in Eq. ([7]) in the rotating frame, this becomes in the RWA 

H 2 = 92 (|tt -| 2 + l a +| 2 ) b^b + -g 2 a-a+b 2 + -g 2 a+a-(r ) 2 = G 2 b^b + -e m b 2 + -(£^) 2 > 


(17) 


(18) 


where 


Gf = VF2 <x± = JgiG±/g 0 
G 2 = (Gp 2 + (G 2 ) 2 
£m = G ! 2 'G ! 2 = g 2 G-G + /gl 


(19) 
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The strength of the parametric modulation is hence given by the quantity e m . In our scheme, ^2/(27r ) 2 - 2 x 1 (T 5 
Hz 2 . We obtain that at the highest pump powers used in the experiment, £ m /27r ~ 300 Hz which is similar to the 
mechanical damping rate, and hence looks like we are close to a parametric instability. However, we find that the 
sideband cooling interaction counteracts parametric modulation, and that the mechanical parametric modulation has 
a vanishing effect for either the quadrature temperatures or output spectra for £ m /27r < 10 kHz. 

Apart from the second-order coupling in Eq. dzt parametric modulation of the mechanical frequency could arise 
from thermal effects [23], with a substantially larger amplitude. We cannot qualitatively analyze this because we do 
not know the thermal picture well, however, we find numerically that the output quadrature spectra would be affected 
in a way which is inconsistent with the measurement (see section [| . The peaks in the spectra would sharpen up from 
the ideal case and would sometimes show a multi-peak structure. Hence, parametric modulation of the mechanical 
frequency in the experiment is rather well excluded. 


Cavity 


In the ideal model together with Eq. parametric modulation of the cavity, or fields rotating at 2 cj c , do not 

arise. However, nonlinearities in the cavity, such as a Kerr effect, uj c => uj c (l + ^KA^A), would induce fields at 2cj c . 
In the rotating frame we obtain 


H c2 = ^K(A*fA 2 = ^Ka-a+{a ] f + ^Ka*_a* + a 2 + K (|a + | 2 + |a_| 2 ) a f a = 
= \ £ c a2 + 7} £ cC ) 2 + K'a’a 


( 20 ) 


Here we also defined the complex modulation amplitude e c . Because the field amplitudes (Eq. 0 ) can be taken as 
real, e c is real as well in case of the Kerr effect. However, as shown in section [J the calibration and squeezing spectra 
indicate a purely imaginary e c (phase shift — 7 r/ 2 ) (section [j). This type of phase shift is indicative of modulation of 
dissipation. 

A possible explanation is nonlinear damping. The extra nonlinear damping term modeled as nonlinear coupling to 
a bath of oscillators [36] in the equation of motion for a X\ equals — rjXi(t) 2 Xi with a real-valued proportionality 
constant 77 . For the equation of motion for the annihilation operator in the lab frame, this corresponds to the term 
\ri(A 2 +2 A^A + (A^) 2 )(A^ — A) added to the rhs. In the present pump scheme, this term becomes — rja-a+a^ which 
corresponds to the required a purely imaginary £ c in the equations of motion, Eq. (23). As typical of microfabricated 
cavities, we indeed observe power-dependent dissipation, but we have not been able to sort out the contributions due 
to different mechanisms, such as two-level system damping and breaking of superconductivity at high currents. This 
would require the understanding, for example, of the time scales of the different processes. 


Input-output modeling 

The model for the transmission setup is depicted in Fig. [5] The following modeling is basically standard input- 
output theory. The asymmetric cavity has the input side where the pumps are applied, with the corresponding 
input coupling rate kei- Energy inside the cavity can dissipate in three separate channels. One is through the input 
coupler. Another channel is the cavity internal losses at the rate ft/, and the third is the output coupler with the 
rate ke 0 ft#*. The latter property enables us to catch nearly twice the signal as compared to a symmetric cavity, 
important for the integration time needed in the measurements. 

The input fields to the cavity are the standard input-output theory quantum fields described by the operators 
&m,o, &in,i for the input and output couplers, and internal losses, respectively. They have the frequency-domain 
correlators (n = i, o, I) 


( a ln,n( LU ) a ^^( ^ 0 ) — n n(u c ) 

= (1 + n n{^c)) 5 


( 21 ) 


other correlators being zero. 

The three baths in general have unequal temperatures T n which give the individual Bose factors 


n n {u c ) = (exp (hw c /k B T n ) - 1) 1 


(22) 
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The internal losses of the mechanics are due to the field bi n which has similar properties as in Eq. (21), Eq. (22), 
with the replacements uj c =>* u; m , n n =>• n^. 

We suppose the i and o baths are at zero temperature which is standard assumption. We also have carefully 
taken care to isolate excess noise from both these ports. The fact that the thermal calibration and sideband cooling 
calibration work well without a dip arising from such noise supports this assumption. The cavity internal bath, 
however, rises up to n/ ^ 0.9 at the highest pump powers. 

The equations of motion, including for completeness the parametric modulations to both the mechanics and cavity, 
as well as pump detunings, are 

a = -iA s a - - ie c o) - i ( G-b + G + 6 f ) + y/nEla in ^(t) + y/KEoQ>in,o{t) + 


= iA s a j - + is* c a + i {G*_b ] + G*+b) + v^ a L,i(*) + V^ a L,o(^ + v^ a L,/W 


(23) 


We solve the frequency-domain equations resulting from Eq. ( [23] ) numerically up to all orders. The response of 
each operator is written as arising from all the input fields, given as an example for the cavity operator as: 


b = -i(A A + G 2 )b ~^b- is m tf - iG+a) - iG*_a + ^b in 
$ = 4 + G 2 )b t — -^-6^ + is^b + iG* + a + + y/'y m b\ n 


a(uj) = m(uj)a in ^(bj) + l(<jj)a\ ni {io) + m 0 (uj)a in ^{uj) + J 0 (w)aT no (w) + 
+ Tni{u)a in j(w) + li(u)a\ n j(u) + q(u)b in (u) + r(u)b\ n (u). 


t 


(24) 


One of the aims is to calculate the quadrature spectra of the mechanical resonator, in particular, the fluctuation 
energy as a function of the quadrature angle O. The spectra of the mechanics is given as 


S x ,&(w) = ((b\oj) exp(iO) + b(u) exp(—iO))(^(— cj) exp(i &) + b(—uj) exp(— iO))) = 

= (bf(co)bf(-co)) exp( 2 iO ) + (b^(cj)b(-cj)) + (b(cj)b^(-cj)) + (b(cj)b(- cj)) exp(— 2 iO). 


(25) 


The correlators are calculated using Eq. ( [24] ) , Eq. ( [21] ) . 

The quadrature variance, directly giving the energy of the quadrature in units of the zero-point energy is given as 


1 f°° 

AX| = — / 2 S x , e (ij)dLj. 

In these units, the total energy of the mode is 

)• 


= 4 ( A ^ 


AX | +7r/2 


(26) 


(27) 


For example, in the ground state this gives n m = 

We use the standard nomenclature that the minimum of AJ@ with respect to O is called the X\ quadrature, 
and the maximum is The corresponding variances are labeled with the shorthand notations and AXf, 

respectively. 

The final measured quantity is the output field of the field leaking out from the output coupler of the cavity, given 
as 


&out,o — \JftEo& { 


(28) 


which can be expressed as a sum of inputs the same way as in Eq. (24). For the first two terms, for example, this 
reads a 0 Ut ^ 0 (uj) = rn 0 Ut: 0 (uj)ai n ^{uj) + l out , 0 {oS)a\ n ^(u;). The spectral density of the corresponding voltage fluctuations 
is labeled Sy(u). This is the quantity presented in the figures, given in units of quanta. The quadrature spectrum 


Sv,e{u) of the voltage fluctuations is given as in Eq. (25) by replacing the mechanics operators with those of the output 
field (Eq. ([28])) . For notational convenience, we use for the spectra of the minimum energy and maximum energy 
quadratures the notations Si (w) and 62 (<^), respectively. The final deductions regarding squeezing in the mechanical 
resonator are based on direct comparison of the measured phase-sensitive output spectrum to the predictions resulting 
from the presented formalism. In principle it is possible to obtain analytical results, but they are too complicated to 
be written down, but nevertheless, allow for numerical analysis. 
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Since theoretical expectations arise from nontrivial analysis, it is integral to make sure the calculations are sound. 
We verified that the results of the theoretical formalism are similar to earlier work in cases where comparison is 
possible. More specifically, the amount of squeezing our formalism gives is the same as in Ref. 15, as well as the total 
output spectrum reduces to their result in the limit of small internal cavity losses. If considering only the regular 
sideband cooling case, analytical results are available in Ref. 28, and our analysis agrees with those results. 

For certain calibrations (section [| we use the direct microwave transmission measurement through the cavity. The 
usual microwave transmission coefficient is given by 

l&iMI = fL 0Ut,0 (^) ffi ^out,o M. (29) 

CALIBRATIONS 

The calibrated parameters of the experiment are listed in Table |IJ In main text it was described how and the 
Bogoliubov mode ratio 1 / tanhr = G_/G + were obtained. In this section, we discuss further details, and how the 
remaining parameters are calibrated. 


TABLE I: List of calibrated parameters. 


symbol 

gI 

Q 

1/tanh r 

7 m 
KEi 
KEo 
El 

As 

A a 


explanation 

effective optomechanical coupling due to the red-detuned pump 

gain of the system (amplifiers plus cabling) following the output port of the cavity 

Bogoliubov mode ratio 

intrinsic damping rate of the mechanical resonator 
external damping rate of the cavity through the input port 
external damping rate of the cavity through the output port 
internal damping rate of the cavity 
symmetric detuning of the pump frequencies 
asymmetric detuning of the pump frequencies 


Effective coupling, system gain: G-, Q 


Let us discuss the basis of the thermal calibration as shown in Fig. 2a in main text. The output spectrum Sv(w) 
divided by system gain Q , when the pump is applied at the red sideband, supposing <C ft, is 


Sv _ ^Ei^l n KEi 7opt7eff 
Q K, 2 ni K (w - UJ m ) 2 + 7 e 2 ff /4 


(<- 2 T nI ) ■ 


(30) 


which is a Lorentzian centered at the cavity frequency. Here, the sideband cooling enhances damping by the amount 

4 G 2 , . 

(31) 


7opt - 5 

ft 

so that the total (” effective”) damping of mechanics is 

7eff 7opt T 7m • 

Here, the effective coupling depends on the photon number induced by the pump: 

G- = goJnZ . 


(32) 


(33) 


We now suppose 7 op t 7 m such that cavity back-action is negligible. Also ni n^, such that we can neglect n/ 
in the second term in Eq. (30). The first term involving ni would cause a shift in the base level. In any case, n/ is 
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essentially zero at the low pump powers. We fit a Lorentzian to each peak, obtaining the amplitude and linewidth for 
each temperature. 

From Eq. (30) we obtain a prediction for the peak area 


Apeak = 2irg^4G 2 _nl = 2^g~^AG 2 _k B T/hw m = 1CT 


(34) 


In cryogenic microwave measurements, Q is poorly known and difficult to measure because of uncertainties in the 
attenuations of cables. For the same reason, the value of G_ at a given input power to the cryostat is difficult to tell. 
In this calibration, therefore, we simply combine all (thus far) badly known parameters into a single coefficient JC. 
We fit Apeak to a linear temperature dependence, obtaining Fig. 2a in the main text. Finally the equilibrium phonon 
number is obtained as = A pea k&B/(fio; m /C) for each data point. 


For calibrating G_ versus generator power, we study the peak width as a function of power (Eq. ([32])) . According 
to Eq. (31), Eq. (32), Eq. (33), the peak width is proportional to the cavity photon number which depends linearly 

directly gives y e ff 


on power. One has to bear in mind that the width of the peak in the output spectrum, Eq. 
only outside the strong coupling regime. The result is shown in Fig. [8j We took into account the effect of around 6 
% change of k versus power. 


xio 4 



generator power (arb) 


FIG. 8: Calibration of effective coupling. The largest power shown corresponds to an effective coupling G m / 2tt ~ 79 kHz. 


For calibrating t/, we use the basic sideband cooling data where example curves shown in black in Fig. 2b-e in main 
text. Each curve (about 20 curves at different power) are simultaneously fitted to theory with the same Q. The theory 
is basically a more accurate version of Eq. ( |30| ) , but allowing for the strong-coupling regime (peaks non-Lorentzian). 
We use the G_ and n ^ just calibrated. The latter means that we know that at low pump powers, the mechanics has 
a given bath temperature. Hence in the fit, we enforce this, but otherwise allow for free and nj (note that these 
n^ and nj are not those of the final squeezing experiment). 


Mechanical linewidth: 7 m 

In Fig. [9] we show how the baths, to which either the cavity and the mechanics are coupled, heat up as a function 
the pump power(s). This heating ultimately limits the squeezing. The bath temperatures are not totally monotonous 
with respect to increasing pump power. We also sometimes observe the bath temperatures to slightly fluctuate during 
a cooldown, presumably due to changes in the configuration of microscopic defects in the sample, varying the thermal 
coupling from the pump into the bath. 

Because the mechanical bath heats up appreciably at the highest pump powers, studying the mechanical linewidth 
at the base temperature of the cryostat (at which temperature the measurements are done) may not be sufficient. 
In Fig. [lO] we show the dependence of the mechanical linewidth as a function of the refrigerator temperature. As seen 
in the figure, the mechanical damping does not depend on temperature more than the scatter of the data points in 
the region < 60 mK. This temperature range covers the bath temperatures at all pump powers. Hence we can use 
the same value for over the whole pump power range. 
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FIG. 9: Bath temperatures as a function of effective coupling. The data points are from the data shown in the main text 
in Fig. 2 (asterisk) and Fig. 3 (circles). As elsewhere, green refers to the BG mode, and black is the sideband cooling, (a), 
mechanics, (b), cavity. 



Temperature (K) 

FIG. 10: Mechanical intrinsic linewidth 7 m vs cryostat temperature. 


Cavity linewidth: keu keo, ki 

The input kei and output ke 0 coupling rates are determined by an electromagnetic simulation in the actual device 
structure. Based on typical spread in the dimensions of fabricated devices, we estimate around 10 % error margins 
for them. The total linewidth k = kei + keo + ki is determined from a fit to the measured transmission in case 
optomechanical effects are not relevant. A subtraction gives the internal decay rate /«/. 


Thermalization of the Bogoliubov mode 


In section D we discussed the thermalization of the bare mechanical mode. Because the BG mode follows the similar 
Hamiltonian, Eq. (13), a similar temperature dependence holds for the energy of the BG mode as well. We calculate 
the green theory curve in Fig. 2a in main text as follows. From the output spectrum based on Eq. ( [25] ) , Eq. (28) we 
numerically calculate predictions for the areas under the sideband peaks separately for both the BG and mechanical 
modes. The green theory curve is by the ratio of these two, a factor of 1.42, above the bare mechanics theory curve. 


Effect of pump detuning: As, A a 

Unless we are operating at the exact sideband co-resonance A± = the resonance condition is lost and physics 

becomes more complicated than simply sideband cooling of the BG mode. The scale at which the ” cavity” detuning 
A s should be correct is k. For the ” mechanics” detuning Aa, the situation is not so clear, but we estimate the scale 
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is between y m and 7 op t, and hence one has to be more cautious. We could fit the detuning to the spectra, but we 


calibrate it independently using the cavity linear response to a probe signal (S21 microwave transmission), see Fig. 11 

It is relatively simple to accurately set the asymmetric detuning, Eq. ( p~3] ) , to zero because we can measure the 
mechanical frequency accurately, and possible changes are expected to be minor, and on the other hand, the cavity 
transmission is sensitive to a non-zero value. Setting the symmetric detuning, however, is more challenging because 
the cavity frequency has some dependence on pump power, and hence we have to first calibrate the cavity frequency 
as a function of pump power. 

An example of the detuning calibration under the double-pump condition is given in Fig. [TlJ For the fit, we ignore 
the parametric effects which have a negligible effect on the transmission. The fit gives As — 11 kHz ±400 Hz, 
A a — 510 Hz ±120 Hz, although the resonance condition was nominally set. Notice that the transmission is quite 
sensitive to both detunings which allows for small error bars. 

At this level of detuning, the affect on squeezing is negligible. For the data in Fig. 2b-e in the main text no careful 
calibration of the symmetric detuning A s was done, but it is estimated to be around 40 kHz based on the shift of the 
BG mode. Even this amount does not cause appreciable degradation of squeezing. 



FIG. 11: Calibration of pump detuning. Plotted is the measured cavity transmission corresponding to Fig. 3d. The theory fit 
is superimposed on the data. 


Parametric modulation 


One can make an independent measurement of the cavity parametric effect by the use of a large detuning; |A±| ^ ft 
(Eq. ®) but such that the center frequency stays at uj c . Expressed with the ”cavity” and ”mechanics” detunings A s 
and A a in Eq. (13), we have A s ~ 0 and A ^4 ft. With the high detuning, optomechanical effects which take place 


within the scale of ft about the sideband resonance, become negligible. 

We carried out such a measurement in the scheme approximately the same as Fig. 3d in main text. We used 
A± ~ uj m ± (27r) • 4 MHz. Although we used the same generator power as in Fig. 3d, the effective couplings are 
smaller because of higher detuning from the cavity; we obtain G-/2tt cs 180 kHz. 

The theory curves in Fig. 12 are calculated for a parametrically modulated oscillator (no optomechanics). We 
obtain a good agreement by using the values \e c \/2it = 42 kHz and the phase arg (e c ) = — 7 t/ 2 . The magnitude is in 
reasonable agreement to what we have for the squeezing spectra with a similar G_ in Fig. 3c. With different values of 
arg (e c ), the pattern would shift by a corresponding amount with respect to LO phase. In particular, the Kerr effect 
of the cavity corresponds to arg (e c ) = 0 which is excluded. 


Different types of parametric modulation 


Here we discuss how various types of parametric modulations would show up in the quadrature spectra. As seen 
in Fig. 13, changes in phase or amplitude, or modulation to the mechanics instead of cavity, each cause characteristic 
features in the spectra, each of which are inconsistent with data. As seen in (a) and (b), the experiment and theory 
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FIG. 12: Calibration of parametric effects, (a) Curves plotted at different LO phase values at 7r/10 steps from 0 to 7r/2, from 
bottom to top. The theory fit is superimposed on the data, (b) The same data plotted as surface, (c) Corresponding theory 
plot. 


clearly show a pattern symmetric with respect to the LO phase 7 t/2 . In particular, the case Fig. 13: corresponds to 


the cavity Kerr effect, causing very asymmetric pattern. In (d), the modulation of the mechanical frequency tends to 
sharpen the peaks. In either case, the effect on squeezing is minor. 


a b 



1.25 1.3 1.35 1.4 1.25 1.3 1.35 1.4 

frequency (Hz) x io 7 frequency (Hz) x io 7 


FIG. 13: Squeezing spectra with different parametric effects. The setup is as in Fig. 3d in main text, (a) Experiment, (b), 
Theory as in Fig. 3d. (c) Same as (b), but with arg (k c ) = 0. (c), Parametric modulation to mechanics: £ c = 0, = 4.2 kHz. 


As mentioned in the main text, it is possible to obtain a good fit of the Bogoliubov mode data (that is, the total 
output spectrum) without imposing any parametric modulation. However, one finds out that far too little phase 
dependence would be predicted for the quadrature spectra. We illustrate this in Fig. 14 The color codes are as in 
Fig. 3d in main text. The parameters are n^ = 91, n/ = 1.07, AXf = 0.79. One has to use an elevated cavity 
internal temperature to explain the enhanced emission of BG mode. In the real case, the parametric modulation acts 
to boost up the emission. The squeezing is only little affected, nevertheless. 


ERROR ANALYSIS 

Next we discuss the error analysis of the final results demonstrating squeezing. We will notice the imprecision is 
dominated by uncertainty of the system gain at the output side. The other significant contributions to the imprecision 
come from those of the mechanics bath thermalization, mechanical intrinsic linewidth, and of the BG mode ratio. 
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FIG 


14: Modeling without parametric modulation. As Fig. 3d in main text, but with £ c = 0. 


A general remark is that the shapes of the BG mode and quadrature spectra are sensitive to most parameters, the 
bath temperatures of both the mechanics and the cavity, and n/ in particular. The third parameter on which 
squeezing significantly depends on and whose independent determination has non-negligible error margins, is the 
mechanical intrinsic linewidth y m . 

We make error estimate for the squeezed quadrature variance using standard error propagation as 


<5 (AX?) = 


N 




(35) 


where the statistical as well as systematic (calibration) uncertainties of variables yi are marked by Syi. The derivatives 


dAXl 

dy 


are obtained numerically. 


The error terms are summarized in Table [TT| If possible, we list both statistical and systematic error for each 
parameter. The quoted statistical errors for the three adjustable parameters nj, £ c are from direct fits to the 
squeezing spectra. For some calibrated parameters, there can be also systematic uncertainties on top of those obtained 
from the fit in question. The statistical errors are quoted as 2 a error bars. The uncertainty for each parameter is 
obtained as the largest of the statistical error and systematic error. 

The uncertainty Sn ^ of the starting temperature of the mechanics bath is particularly important to consider. At 
low pump powers, an error in the bath temperature calibration would show up as a scaling factor in the Lorentzian 
peak, easily going unnoticed. The real situation, however, is better because towards the strong-coupling regime, the 
starting temperature has a more complicated influence than just scaling the peaks. The linear fit to Eq. ( [34] ) (plotted 
in Fig. 2a) gives a negligible uncertainty, only ^ 1 %. A systematic error estimate is based on the fact that the data 
points at T <20 mK, which region marks the base temperature, show more scatter than those at higher temperature, 
probably marking fluctuations in thermalization. We get another estimate for dn^ from the scatter. The standard 
error of the mean gives Sn ^ ~ 3, or approximately 10 % uncertainty. One more way to obtain Sn ^ is from the direct 
statistical error of the fits to the quadratures, yielding Sn ^ ~ 1.2. 

The uncertainty of the fitted value of in the data shown in Fig. le in main text is d^ m /27r ~ 14 Hz. This 
value is considered as the systematic uncertainty. Because is not an adjustable parameter in the end (when fitting 
to the quadrature spectra), it does not have a statistical uncertainty because of the way we carry out the analysis. 
Nevertheless, one might ask the question of what if would have a statistical uncertainty. For the latter purpose, 
we can make a test where in the fits we in fact consider it as an adjustable parameter. Some of the true adjustable 
parameters were fixed during this test for numerical stability. The direct statistical error this way gives an estimate 
^7 m /27r ^ 5 Hz, which is negligible in comparison to the systematic uncertainty. 

The uncertainties of n/ and e c are obtained as a direct statistical error of the fits to the quadratures. 

The uncertainty of the cavity linewidth k is obtained from Lorentzian fits to the measured cavity transmission. The 
statistical error from a given fit is negligible, but we obtain some scatter between different repetitions over a longer 
time period. Standard deviation of the latter gives 5k/2tt ~ 10 kHz. Based on typical spread in the dimensions of 
fabricated devices, we estimate around 10 % error margins for kei and keo, and hence that for ki is the sum of all 
three. We note the squeezing is very insensitive to how n is distributed between different channels. 


The confidence level for G_ when fitted to Eq. (32) is ±2 %. 
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TABLE II: Error budget. The numbers are for the optimum squeezing, i.e., Fig. 3d in main text, G_ = 235 kHz. The last 
column tells the contribution of the error term to the total uncertainty of the squeezing. The frequencies are given in linear 
frequencies (Hz). 


parameter yi 

value 

statistical 8yi 

systematic 8yi 

statistical 5yi d Q* 1 

systematic Syi d g^ 1 

percentage 

T 

91 

1.2 

3 

0.006 

0.015 

11 

ni 

0.81 

0.02 

- 

0.005 

- 

4 

£c 

56 kHz 

1.2 kHz 

- 

0.001 

- 

1 

G- 

235 

- 

5 kHz 

- 

0.01 

7 

G 

arb 

2 % 

11 % 

0.01 

0.06 

41 

1/tanh r 

1.43 

- 

0.1 

- 

0.014 

10 

7m 

330 Hz 

5 Hz 

14 Hz 

0.006 

0.016 

11 

KEi 

50 kHz 

- 

5 kHz 

- 

10 -3 

2 

KEo 

270 kHz 

- 

30 kHz 

- 

0.0039 

3 

Kl 

330 kHz 

- 

40 kHz 

- 

0.016 

11 

A s 

11 kHz 

400 Hz 

- 

10 -5 

- 

0 

A a 

510 Hz 

120 Hz 

- 

10- 5 

- 

0 




For Q, the direct statistical error is small, but we also have to propagate errors from the parameters G_ and n ^ 
which are used to obtain Q. The error due to G_ is negligible, but that due to equals the confidence level of n^, 
around 10 %. In order to obtain the corresponding error in squeezing, we enforce a change in Q from the calibrated 
value and re-run the fits. As seen in Table [TT| the gain error dominates the squeezing uncertainty. 

Uncertainty for the pump power ratio 1/ tanhr = G_/G + is obtained as systematic uncertainty of the amount 0.1 
equaling the drift mentioned in the main text. 

Combining the terms in Table [IT] we obtain 


S(AX^) -0.08, 


(36) 


which is the final uncertainty estimate. 

For smaller pump powers, the uncertainties are larger because the peak in the cold quadrature is lower and signal- 
to-noise hence smaller and statistical uncertainties begin to dominate. 

We can also demonstrate that the theory curves are rather sensitive to the parameters. In Fig. 15 we show plots as 
in Fig. 3 d in main text, but one of the important parameters is intentionally offset by 20 % from the nominal value 
to the direction which degrades squeezing, while other parameters are kept unchanged. In each case, the agreement 
is clearly worsened, while squeezing is affected around 10 %. 



frequency (Hz) X 10 7 




frequency (Hz)' x10 7 ' frequency (Hz) x -| 0 7 


FIG. 15: Tests of sensitivity to parameters, (a), = 109, (b), ni — 0.97, (c), 7 m /27r = 396 Hz. 

















